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ABSTRACT 

In order to constrain and possibly detect unusual physics during inflation, we 
allow the power spectrum of primordial matter density fluctuations, P in (k), to 
be an arbitrary function in the estimation of cosmological parameters from data. 
The multi-resolution and good localization properties of orthogonal wavelets 
make them suitable for detecting features in P in (k). We expand Pi n {k) directly 
in wavelet basis functions. The likelihood of the data is thus a function of the 
wavelet coefficients of Pi n (k), as well as the Hubble constant Hq, baryon density 
Qt,h 2 , cold dark matter density Q c h 2 and the reionization optical depth 
flat ACDM cosmology. We derive constraints on these parameters from Cosmic 
Microwave Background (CMB) anisotropy data (WMAP, CBI, and ACBAR) and 
large scale structure (LSS) data (2dFGRS and PSCZ) using the Markov Chain 
Monte Carlo (MCMC) technique. The direct wavelet expansion method is dif- 
ferent and complimentary to the wavelet band power method of Mukherjee & 
Wang (2003a,b), and results from the two methods are consistent. In addition, 
as we demonstrate, the direct wavelet expansion method has the advantage that 
once the wavelet coefficients have been constrained, the reconstruction of Pi n (k) 
can be effectively denoised, i.e., Pi n (k) can be reconstructed using only the coef- 
ficients that, say, deviate from zero at greater than la. In doing so the essential 
properties of Pi n {k) are retained. The reconstruction also suffers much less from 
the correlated errors of binning methods. The shape of the primordial power 
spectrum, as reconstructed in detail here, reveals an interesting new feature at 
0.001 < k/Mpc' 1 < 0.005. It will be interesting to see if this feature is con- 
firmed by future data. The reconstructed and denoised Pi n (k) is favored over the 
scale- invariant and power-law forms at > la. 



1. Introduction 

Inspite of the success of the inflationary paradigm (Guth 1981; Albrecht & Steinhardt 
1982; Gott 1982; Linde 1983; Kolb & Turner 1990; Hu 2003; Peebles & Ratra 2003), we are 
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still in search of the correct model of inflation. The simplest models of inflation predict a 
power-law primordial matter power spectrum (for example, Linde (1983); Freese, Frieman, 
& Olinto (1990); La & Steinhardt (1991)). However, there are many other viable models of 
inflation which predict primordial power spectra which cannot be parametrized by a simple 
power law (for example, Holman et al. (1991ab); Wang (1994); Randall, Soljacic, & Guth 
(1996); Adams, Ross, & Sarkar (1997); Lesgourgues, Polarski, & Starobinsky (1997)). These 
can represent unusual physics in the very early universe (for example, see Chung et al. 
(2000); Enqvist & Kurki-Suonio (2000); Lyth, Ungarelli, & Wands (2002); Feng & Zhang 
(2003)), and then assuming a power-law primordial matter power spectrum could lead to our 
missing the discovery of the possible features in it and erroneous estimates of the cosmological 
parameters (Kinney 2001). 

With results from WMAP, CMB data continue to be consistent with inflation (Ben- 
nett et al. 2003, Spergel et al. 2003), having even detected the predicted anti-correlation 
between CMB temperature and polarization fluctuations near I of 150. While the derived 
cosmological parameter constraints are consistent with previous constraints from indepen- 
dent and complementary observations, Peiris et al. (2003) report that the data, both with 
and without the use of complementary LSS data, seem to indicate a preferred scale in the 
primordial power spectrum. They conclude that the data suggest at the 2a level, but do not 
require, that the scalar spectral index runs from n s > 1 on large scales to n s < 1 on small 
scales. Their analysis was done assuming a primordial power spectrum of power-law form 
about a pivot scale of 0.05 Mpc" 1 (0.002 Mpc -1 ) when using scalar (both scalar and tensor) 
modes. 

Now that such precision measurements of CMB anisotropy have been made, we aim 
to allow the primordial power spectrum to be an arbitrary function, independent of any 
particular inflationary model, and use the data to place constraints on its shape. WMAP 
data underscore the importance of such measurements suggested earlier by Wang, Spergel, & 
Strauss (1999); Wang & Mathews (2002) who used conventional top-hat binning and linear 
interpolation respectively in the parametrization of Pi n (k). Mukherjee & Wang (2003a) 
parametrized Pi n (k) in terms of wavelet band powers, and arrived at a detailed measurement 
of Pinik) using pre- WMAP CMB data. A preliminary indication of a feature was reported. 
The method was shown to work better than the linear-interpolation method of Wang & 
Mathews (2002). Mukherjee & Wang (2003b) (hereafter MW03b) present constraints on 
Pi n (k) from WMAP data using the wavelet band powers method, and also the tophat binning 
method for comparison. It was found that the reconstructed spectrum deviates from scale- 
invariant as well as power-law models by approximately la. While the data do not rule out 
simple parametrizations such as scale- invariance, as also concluded by Bridle et al. (2003), 
Barger, Lee & Marfatia (2003), Seljak, McDonald & Makarov (2003), there is indication of 
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excess power or a preferred scale at k ~ 0.01 Mpc 1 (Peiris et al. 2003, Spergel et al. 2003, 
MW03b). Wang & Mathews (2002); Mukherjee & Wang (2003a) show that pre-MAP data 
indicated a similar feature. 

It is worth investigating this issue further. Since data in CMB multipole space is essen- 
tially being used to constrain P in (k), its band powers as estimated from the data become 
correlated with one another (because of the cosmological model dependent non-linear map- 
ping between wavenumber k and multipole /, see Tegmark & Zaldarriaga 2002). The Pi n (k) 
parameters in any model parametrization are correlated with the cosmological model pa- 
rameters (that describe the evolution of structure from the initial conditions to the present). 
These correlations can make it tricky to efficiently constrain Pi n (k). 

In this paper, we explore a more direct method of using wavelets to probe Pi n (k); we 
parametrize Pi n (k) by the coefficients of its wavelet expansion. We obtain constraints on 
these wavelet coefficients and cosmological parameters from current CMB and LSS data. 
This direct wavelet expansion method is different and complimentary to the wavelet band 
power method of MW03a,b. 2 It has the advantage that once the wavelet coefficients have 
been constrained, the reconstruction of Pi n {k) can be effectively denoised, i.e., Pi n {k) can 
be reconstructed using only the coefficients that deviate from zero at greater than la. In 
doing so the essential properties of Pi n {k) are retained. Because wavelets are adaptive, and 
the wavelet basis complete, any large scale as well as localized features in Pi n (k) can be 
efficiently picked up by a small number of coefficients. Such properties make wavelets useful 
tools in signal reconstruction. This method is thus used to reconstruct Pi n (k) as an arbitrary 
function over the entire range in k that cosmological data are sensitive to, and since it can 
be reliably reconstructed from a small number of wavelet coefficients, it suffers much less 
from the correlated errors of binning methods. 

We describe our method in Sec. 2. We present our results in Sec. 3. Sec. 4 contains a 
summary and discussion. 



2 The direct wavelet expansion method uses a set of orthogonal basis functions that are very different from 
those used in binning methods; the resultant constraints on Pi n (k) can be thus checked for consistency. This 
can be significant because parameter correlations, which are important for this problem, are very different 
for the different sets of basis functions. We found the parameter correlations to be the smallest in the direct 
wavelet expansion method. 
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2. 



Method 



We expand the primordial power spectrum Pi n (k) 3 directly in wavelet basis functions, 



j=0 1=0 

where ipjj are the wavelet basis functions, constructed from the dilations and translations of 
a mother function tp(k), via (Daubechies 1992; Press et al. 1994) 



The resulting wavelet bases are discrete, compactly supported and orthogonal with respect 
to both the scale j and the position / indices, and hence their coefficients provide a complete 
and non-redundant (hence invertible) representation of the function P in {k). The scale index 
j increases from to J — 1, and wavelet coefficients with increasing j represent structure in 
Pin(k) on increasingly smaller scales, with each scale a factor of 2 finer than the previous 
one. The index /, which runs from to 2 J — 1 for each j, denotes the position of the 
wavelet basis ipjj within the jth scale. (Note that the properties of the basis functions are 
completely different from the basis functions in binning methods). Therefore the wavelet 
coefficients in Eq.(l), together with the scaling function coefficient which represents the 
mean of the function, (denoted a 0j o), total to 2 J in number; correspondingly the function 
Pin(k) is reconstructed at 2 J points. We have considered the case when these sample points 
ki, over which the wavelet basis is defined and at which the function Pi n {k) is reconstructed, 
are equally spaced in log(k). We use 16 wavelet coefficients to represent Pi n (k) in the 
range 0.0002 < A;/(Mpc _1 ) < 0.2 that the data are sensitive to, with the spectrum for 



k < 0.0002 Mpc" 1 set to P in (k = 0.0002 Mpc" 1 ) and the spectrum for k > 0.2 Mpc" 1 set to 
P m (A; = 0.2Mpc- 1 ). 



Thus, when the sample points are equally spaced, the coefficient bjj measures the signal 
in Pi n (k) on scale L/2- 7 , and centered at k = /L/2 J , where L is the full length over which Pi n (k) 
is being expanded. The simultaneous and adaptive localization of the wavelet basis functions 
in position and frequency space make them useful tools for analyzing a variety of data. The 
information content of most signals can be efficiently encoded in wavelet space; hence the 
wavelet representation of a signal is often sparse and allows for significant compression (as 
coefficients below a threshold can be thrown away with very little loss of information making 



3 For a power-law primordial power spectrum, P n = 2.95 x 10 9 A (k/k ) ns 1 (Spergel et al. 2003). For 
an arbitrary primordial power spectrum, we define Pi„(k) = 2 95^10-0 ■ 
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them important tools in denoising). For example a scale-invariant Pi n (k) requires only two 
non-zero coefficients (the scaling function coefficient, and the wavelet coefficient of scale 
0), and when attempting to reconstruct this function from noisy data all other wavelet 
coefficients would be zero within error bars. If some other wavelet coefficients are non-zero 
then that implies additional structure in the function. 

Thus our goal is to constrain Pi n {k) by estimating its wavelet coefficients, together 
with cosmological parameters, from a likelihood analysis of cosmological data. We do not 
impose any particular inflation model dependent parametrization on P in (k), and allow it to 
be an arbitrary function. After constraining the wavelet coefficients, marginalizing over all 
other parameters considered, we can choose to reconstruct Pi n {k) using only those wavelet 
coefficients that are significantly constrained. In doing so all the important features of Pi n (k) 
are retained and one is essentially denoising. Thus in this construction features in any Pi n (k) 
can be localized in scale j and position I with the help of a few coefficients, and the resulting 
reconstruction also suffers much less from the correlated errors of binning methods. 

This direct wavelet expansion method is very different in nature from the wavelet band 
power method used in MW03a. In the wavelet band power method of MW03a, the primor- 
dial power spectrum is parametrized in terms of wavelet band powers; the window functions 
relating the wavelet band powers to the primordial power spectrum are the squares of the 
Fourier transforms of the wavelet basis functions ( apart from a j dependent normalization 
factor). Wavelets that are more smooth and hence less localized in pixel space and cor- 
respondingly more localized in Fourier space, such as the Daubechies 20 wavelet are the 
suitable wavelets to use for this parametrization. In the direct wavelet expansion method we 
use the coefficients of a wavelet expansion of Pi n {k) to characterize it; the window functions 
in this case are the mutually orthogonal wavelet basis functions iftjj. Thus here we need to 
use wavelets that are most localized in k space, which is the pixel space here, such as the 
Haar and the Daubechies 4 (hereafter D4) wavelet. Among such wavelets different wavelets 
give very similar results. Other differences between the two methods are that in the wavelet 
band power method the bin locations and their separations are not arbitrary, as explained 
further in MW03a (see also footnote 7 here), whereas here we can probe the function down to 
arbitrary resolution by reconstructing more coefficients. The correlations among the recon- 
structed parameters are different and smaller in the direct expansion method as compared 
to those in binning methods, and the reconstruction here can be effectively denoised; these 
properties have been elaborated on earlier in this section. 

We also make use of the mutual independence of the wavelet basis functions in computing 
the projections of the CMB angular power spectrum Cz's and the matter power spectrum 
P{k) for each wavelet basis, thus avoiding calling CMBFAST (Seljak & Zaldarriaga 1996) 
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or CAMB (Lewis, Challinor, & Lasenby 2000), when only the wavelet coefficients are varied 
for the same cosmological parameters. This leads to a significant computational speed up. 
We write 4 

-P m {k) |A Ti (A;,r = r )| 2 
= / T ^M\^T l (k,T = T )\ 2 

3,1 

= E^^W' ( 3 ) 

where the cosmological model dependent transfer function Axi(k, t = r ) is an integral over 
conformal time r of the sources which generate CMB temperature fluctuations, r being the 
conformal time today, and s represents cosmological parameters other than the bj/s. The 
C/''(s)'s are the wavelet projections of the Cj's, computed here using CAMB 5 . 

Similarly the wavelet projections of the evolved matter power spectrum P(k) are com- 
puted for LSS data as follows, 

P(k,{b jtl },s) = P m (k)*T(k) 2 

= E^(*)*n*) J 

ee J2hiP j ' l (k,s) (4) 

3,1 

where T(k) is the matter transfer function. 

Since the theoretical predictions to be compared with CMB and LSS data both depend 
on P in (k), and the cosmological parameters, they can together be used to place robust 
constraints on cosmological models and their initial conditions. Note however that as pointed 
out by Elgaroy, Gramann, & Lahav (2002), it is hard to detect features in the primordial 
power spectrum at k < 0.03/iMpc -1 using LSS data only. 



4 In computing these wavelet projections, we note that tpk,i are defined only at a finite number of points 
(here 16) and its value at any intermediate point is obtained by linear interpolation. This is appropriate; 
since we are using a discrete transform and correspondingly estimating a function at a finite number of 
points, no additional information at smaller scales is either assumed or included. We have checked that no 
interpolation errors are introduced at this stage. 

5 http://camb.info/ 
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We use the Markov Chain Monte Carlo (MCMC) technique, illustrated for example 
in Lewis & Bridle (2002); Verde et al. (2003) to estimate the likelihood functions of the 
parameters. This technique is necessitated by the large number of parameters being varied 
here, and at the same time is free of the interpolation errors expected in the much slower 
grid based methods. At its best, the MCMC method scales approximately linearly with 
the number of parameters. See Neil (1993) for a review, and Hannestad (2000); Knox et 
al. (2001); Rubino-Martin et al. (2002) for other applications of this method to CMB data 
analysis. An MCMC sampler based on the Metropolis-Hastings algorithm has been made 
available in the software package CosmoMC (Lewis & Bridle 2002). 

The MCMC method is used to trace the full posterior distribution of the parameters in 
a flat ACDM cosmology assuming wide uniform priors of 0.4 < h < 1.0, where h = H /(100 
kms -1 Mpc -1 , 0.005 < Vl b h 2 < 0.1 and 0.1 < VL c h 2 < 0.99 on the cosmological parameters. 
The two wavelet coefficients that set the level of Pi n {k) are allowed to vary between 
and 5, while the remaining wavelet coefficients are allowed to vary between ±2.5, with 
the requirement that the resulting Pi n (k) be positive definite 6 . We use a weak prior on 
the age of the universe of to > 10 Gyrs. We have run several chains, and have checked for 
convergence and mixing. From the MCMC samples thus obtained, the marginalized posterior 
distributions and confidence limits of the parameters are estimated. 



3. Results 

We use CMB temperature anisotropy data from WMAP (Bennett et al. 2003), com- 
plimented at / > 800 by data from CBI (Pearson et al. 2002) and ACBAR (Kuo et al. 
2002) upto an l max of 2000, marginalizing analytically over known beamwidth and calibra- 
tion uncertainties (Bridle et al. 2002). We use LSS data from 2dFGRS (Percival et al. 2002; 
the redshift space power spectrum is assumed to be proportional to the real space power 
spectrum on the large scales used) and PSCZ (Hamilton & Tegmark 2002; decorrelated real 
space power spectrum) galaxy redshift surveys over linear scales. A linear bias is assumed 
and analytically marginalized over (Bridle et al. 2002). LSS data are thus effectively being 
used to constrain the shape and not directly the amplitude of the matter power spectrum. 

We assume Gaussian adiabatic scalar perturbations in a flat universe with a cosmological 
constant. We do not use tensor modes in this paper, since current data are not sensitive to 
tensor contributions. With freedom allowed in the shape of the primordial power spectrum, 



6 Note that for a scale-invariant Pi n (k) of unity the two coefficients that set the level are each 2.828 while 
the remaining 14 coefficients are 0. 
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the temperature anisotropy data themselves cannot distinguish well between different values 
of the reionization optical depth r Ti . Hence we simply make use of the WMAP constraint of 
r ri = 0.17 ± 0.04, derived for a flat ACDM cosmology from WMAP's TE polarization data 
(Kogut et al. 2003), where the error estimate includes systematic and foreground uncertain- 
ties. We marginalize the distributions of the other parameters over a Gaussian distribution 
in r ri centered at 0.17 with a standard deviation of 0.04. For this, we calculate the likelihood 
distributions of the wavelet coefficients as well as the Hubble constant H , baryon density 
fib h 2 and the cold dark matter density Q c h 2 at 5 select values of r ri and approximate the 
marginalization over r r j using the 5-point Gauss-Hermite quadrature summation. 

Fig.l shows the marginalized Id distributions of the parameters, obtained from an 
MCMC based likelihood analysis of CMB and LSS data, after marginalizing over the as- 
sumed distribution for r ri . All the parameters are well constrained, lying well within the 
imposed priors. 

Fig.2 shows the reconstructed Pm(k) (solid curves) and its la confidence region (dotted 
curves), obtained from the mean values and standard deviations of the wavelet coefficients. 
The uncertainty in P in (k) has been calculated using: 



Pin {&i ) 



= iv 

N ^ 



N L j,l 



dPin\ki) 



(5) 



where the average is over the MCMC samples, bjj denotes the mean value of the wavelet 
coefficient, and dP ^ k ^ = tpjj(ki); correlations between the wavelet coefficients have thus 
been accounted for in computing the uncertainty in P in {k). The upper panel shows the 
reconstructed Pi n (k) using all wavelet coefficients. These constraints are shown overlaid on 
the constraints obtained by MW03b using the wavelet band power method. 7 We see that 
while the constraints from the two methods are completely consistent, a new feature shows 
up at large angular scales around 0.001 < A;/Mpc _1 < 0.005. The lower panel shows the 
spectrum reconstructed using only those coefficients that deviate from at > la. These 
constraints are shown overlaid on the constraints obtained using the tophat binning method 
with bands carefully positioned closer together so as to be able to pick up the feature at large 
angular scales around 0.001 < /c/Mpc -1 < 0.005, and placed further apart in the remaining 
parts of the spectrum. The feature shows up consistently. 



7 Note that in the wavelet banding method the locations and the widths of the bands are not arbitrary. 
As explained in Mukherjee & Wang 2003a, the band power window functions in this case are the modulus 
squared of the Fourier transforms of the wavelet basis functions. By definition then, since each scale is a 
factor of 2 finer than the previous one, the bands lie log w 2 apart in k. 
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The coefficients that are significant at the greater than la level are 60,0, &i,o, &2,o, 
62,1 and & 32 , together with the scaling function coefficient. Notice that with 10 of the 
16 wavelet coefficients dropped, the spectrum reconstructed in panel (b) still has all the 
essential features of the spectrum shown in panel (a). The rise in power that shows up 
near k ~ 0.01 Mpc -1 mimics a power-law model with possible running and is what was 
reported in MW03b and by WMAP. In this paper we report a localized new feature at 
0.001 < k/Mpc' 1 < 0.005 that would be missed by a power-law parametrization, and by 
other model-dependent parametrizations. (We note that section 5 of Peiris et al. (2003) 
seems to indicate the same large scale feature obtained here, but they used the model de- 
pendent parametrization of Adams, Cresswell, & Easther (2001) and constrained the Pi n (k) 
parameters at a fixed set of cosmological parameters.) This new feature would also be missed 
by a model-independent binning method if the bin separation were not small enough at the 
appropriate locations. 8 

In the direct expansion method the constraints obtained simultaneously on the cosmo- 
logical parameters are h = 0.59 ± 0.09, Q b h 2 = 0.018 ± 0.004 and Q c h 2 = 0.143 ± 0.019. In 
all the constraints presented here we have used the parameter mean values together with 
marginalized la standard deviations. 9 

The xtff °f the best fit model in the direct wavelet expansion method is 1032.45 for 
a total of 964 data points, when constraining 16 wavelet coefficients and 3 cosmological 
parameters with an external prior on r ri being imposed and marginalized over. 10 A scale- 
invariant model with 5 constrained parameters has a xtff °f 1044.33 for the same data. One 
can see that this model is equivalent to holding fixed the 14 wavelet coefficients that are 
responsible for any deviations from scale-invariance. Thus the difference in x 2 between the 
two models has a ^-distribution with 14 degrees of freedom. By computing the probability 
of obtaining such a x 2 we are able to compare the two models and find that the scale- 
invariant model is disfavoured at ~ 0.6a. The xlff °f the best fit model in the direct 
wavelet expansion method using only the 6 significantly constrained wavelet coefficients is 
1039.52 for 9 constrained parameters. This is favored over the scale-invariant model at 



8 Note that a uniform binning with a small bin size can lead to ovcrsampling of the data, which would 
result in large correlated error bars for the estimated bin amplitudes, thus burying any possible features 
(Bridle et al. 2003). 

9 As discussed in Lewis & Bridle (2002), the MCMC samples from the posterior do not provide accurate 
estimates of parameter best-fit values, because in higher dimensions the best-fit region typically has a much 
higher likelihood than the mean but occupies a very small fraction of parameter space. 

10 Thc number of data points in the different data sets used are 899 (WMAP), 4 (CBI), 7 (ACBAR), 32 
(2dFGRS) and 22 (PSCZ). 
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> lcr. 11 Note that in our analysis the number of constrained parameters are large and the 
significance of deviations from scale- invariance are small because we allow Pi n (k) to be an 
arbitrary function over the entire k range. While the deviation in the shape of Pi n (k) from 
scale-invariance is not found to be significant at this stage, the primordial power spectrum 
reconstructed here seems to indicate intriguing features that could indicate either unusual 
inflationary physics or hidden systematics in current CMB data. 

4. Conclusions 

It is important to reconstruct the primordial power spectrum as a free function from 
cosmological data, the aim being to test the assumptions that are usually made about Pi n (k), 
and to constrain it model-independently so as to be able to detect any features in it that may 
be signatures of unusual physics during inflation. We have presented a new and powerful 
method to extract the primordial power spectrum Pi n {k) as an arbitrary function from 
observational data — the direct expansion of Pi n {k) in wavelet basis functions. 

The wavelet basis functions afford an efficient (sparse) and non-redundant (invertible) 
expansion of a function, with coefficients that encode information about the structure of 
the function on different scales, and with adaptive position resolution. These properties 
make wavelets suitable tools for directly probing the primordial power spectrum for features. 
Reconstructing Pi n (k) using only the wavelet coefficients that deviate from at > lcr amounts 
to denoising the resultant power spectrum while retaining all its important features (Fig 2). 
The direct wavelet expansion method also yields smaller parameter correlations than the 
binning methods. 

The Pin(k) reconstructed in this paper (Fig.2) is favored over the scale-invariant and 
power-law parametrizations at > la. The results of this paper are completely consistent 
with the spectrum reconstructed in MW03b using the wavelet band power method. These 
two methods are very different in nature and complementary (see section 2). 

The direct wavelet expansion method has enabled us to find indication of a new localized 
feature in Pi n (k) at large scales around 0.001 < /c/Mpc -1 < 0.005. It will be interesting to 
see whether future data will confirm the existence of this feature. 



11 The Xeff °f the best fit power-law model for these data is 1043.68 with 6 constrained parameters. In 
order to compare this model with the wavelet coefficients based model one would need to use a more involved 
method, possibly that of Baycsian evidence. 
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Fig. 1. — The figure shows the marginalized Id distributions of all the parameters used; 
16 wavelet coefficients and 3 cosmological parameters. A Gaussian distribution of r has 
been marginalized over (see text). The mean values and la standard deviations on the 
cosmological parameters are h = 0.59 ±0.09, Vt b h 2 = 0.018 ±0.004 and Vt c h 2 = 0.143 ±0.019. 
Corresponding constraints on the wavelet coefficients imply constraints on Pi n (k) shown in 
Figure 2. 
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Fig. 2. — reconstructed using all wavelet coefficients (solid curve) with its la confidence 
region (dotted curves). These constraints are shown overlaid on the constraints obtained by 
MW03b using the wavelet band power method, (b) Pi n {k) reconstructed using only those 
coefficients that deviate from at > la (solid curve) with the corresponding la confidence 
region (dotted curves). These constraints are shown overlaid on the constraints obtained 
using the tophat binning method with bands carefully positioned closer together so as to 
be able to pick up the feature at large angular scales around 0.001 < fc/Mpc -1 < 0.005, 
and placed further apart in the remaining parts of the spectrum. The dashed line is the 
scale-invariant spectrum that best fits these data. 
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